PERIODICALLY GENERATED PROPAGATING PULSES* 



Abstract. Ceftaib. eI|tiiiMIJiife,A4fi,tMn1Kgfffl)BiiaMSi.inte]S^ pulses 
of a field-like unknown while a current-like unknown oscillates periodically with time. A general 
asymptotic theory of this phenomenon, the generalized Gunn effect, has been found recently. Here 
we extend this theory to the case of nonlinearities having only one stable zero, which is the case for the 
usual Gunn effect in n-GaAs. Our ideas are presented in the context of a simple scalar model where 
the waves can be constructed analytically and explicit expressions for asymptotic approximations 
can be found. 
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1. Introduction. The Gunn effect is the periodic oscillation of the current in a 
passive external circuit attached to a dc-voltage biased semiconductor whose electron 
drift velocity has a single maximum as a function of the electric field (and therefore 
the curve of electron velocity versus electric field has negative slope for field values 
on a certain interval, a fact called negative differential mobility) |22ll25j . During each 
period of the oscillation, a pulse of the electric field is created at the injecting contact, 
moves through the semiconductor, and is annihilated at the receiving contact. While 
originally observed in bulk n-GaAs samples, similar current oscillations, mediated by 
pulse dynamics in dc voltage biased semiconductors, have been found in many mate- 
rials, several of which lack negative differential mobility 1 . Instead, other processes 
(impact ionization at impurities P'4. nonlinear capture coefficients |21) . nonlinear 
recombination processes, etc) may be responsible for a current vs. local electric field 
characteristic curve displaying a local maximum followed by a region of negative slope 
(negative differential conductivity). 

Propagation of pulses occurs in many systems of interest in Biology, Physics, . . . : 
morphogen pulses or spikes in activator-inhibitor reaction-diffusion systems modeling 
cell development or chemical reactors (16i il2l lOl 123) . propagation of nerve impulses 
along myelinated or unmyehnated fibers j^EHlElj pulse propagation through car- 
diac cells USIj calcium release waves in living cells semiconductor superlattices 
[SJ I26| and oscillatory instabilities of the current in bulk semiconductors with an 
N-shaped current-field characteristic [Tl 1221 1^ . These distributed systems can be 
spatially discrete or continuous and can be described by a variety of model equations. 
Sometimes a pulse is created from an appropriate initial condition and it reaches a sta- 
ble shape, moving uniformly until it arrives at a boundary. Sometimes understanding 
pulse dynamics is the key to describing a more complicated evolution of the system. 
A good example is the Gunn effect in bulk semiconductors. 

While Gunn-like instabilities have been known for a long time, only recently have 
pulse annihilation and creation at boundaries been studied by asymptotic methods 
|14l Ej. These theories treated the case in which the relevant nonlinear source 
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Fig. 1.1. Self-sustained oscillations of the 'current' J{t) (see inset, which shows one period of 
J{i)) mediated by pulses of the 'field' u{x,t), which are solutions of the model described in Section 
2. The pulse profiles correspond to the times marked in the inset. Case (a): bistable source. Case 
(b): source with a single stable zero. 



term has two stable zeros. Then there are stable wavefronts joining these two zeros, 
and a moving pulse is a flat region of high held bounded by two wavefronts. The 
pulse changes its size if its leading and trailing wavefronts move at different speeds. 
Fig. II. ll shows the time-periodic oscillation of a 'current' J(t) accompanied by the 
repeated generation and motion of flat-top pulses of an 'electric field' u(x,t), which 
are solutions of model equations described in Section |21 Fig. ILlf a) corresponds to 
the case of a bistable nonlinear source. The asymptotic analysis of the Gunn effect 
is based on the dynamics of such pulses 0. However, the source term (velocity- 
field characteristic curve) in very relevant materials, such as bulk n-GaAs [221 1 wide- 
miniband GaAs/AlAs superlattices |ll[ll5j . and semi- insulating GaAs [21], does not 
have two stable zeros. Instead, the nonlinearity has a single stable zero, so the previous 
theories, based on two moving wavefronts, are not valid. Fig. ILlT b) shows the self- 
sustained oscillations corresponding to this case. Can we find an asymptotic theory of 
pulse mediated oscillations in this case? The answer is yes, as we show in this paper. 

To emphasize that our analysis applies to a class of models, we shall analyze a 
simpler problem than that of the Gunn effect. We will study a nonlinear reaction- 
diffusion-convection equation with an integral constraint 6 and a piecewise linear 
source term. Then the pulses can be constructed analytically, their size can change 
as they move, and explicit expressions for the asymptotic approximations can be 
found. Such a construction was used by Rinzel and Keller for the FitzHugh — Nagumo 
equation [20]. For Kroemer's model, a piecewise linear electron velocity was used to 
calculate the exact form and speed of a steadily propagating pulse ^Hl ■ 

The outline of the paper is as follows. Section |21 presents our simplified model. 
It also reviews the key ideas of previous asymptotic theories, valid for an N-shaped 
nonlinearity with two stable zeros. Section |31 discusses the construction of stationary 
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solutions, and the kinematics of wave fronts, in the hniit of long samples. When 
the nonlinearity is piecewise linear, the pulses can be found analytically. Section 0| 
discusses the dynamics of a single pulse moving from the injecting to the receiving 
boundary. We show that the pulse changes its form and speed adiabatically, following 
the instantaneous value of the current. In Section [5l we complete our description 
of the oscillations by explaining what happens when the pulse reaches the receiving 
boundary, and how a new pulse is created at the injecting boundary. SectionElcontains 
our conclusions. The Appendices are devoted to technical matters. 

2. Simple scalar model. The model consists of a one-dimensional nonlinear 
parabolic equation for u{x, t) (the "electric field") with an unknown forcing term J(t) 
(the "current"). There are also an integral constraint (the "voltage bias condition"), 
as well as boundary and initial conditions: 

(2.1) m+''d-. = d^'-'-'^''^' 

1 --^ 



(2.2) - / u{x,t)dx = (I), 

L Jo 

(2.3) u(0,t)^pJ{t), p>0, _(L,i) = 

(2.4) u{x, 0) = f{x) > 0, Q<x<L. 

Here K, p and </> are positive constants. 

This model becomes the Kroemer model of the Gunn effect if the advection con- 
stant K in Eq. (EH) is replaced hy g{u) H H d 02 0^1 ■ For the Kroemer model, 
existence and uniqueness of solutions have been studied by Liang j^. Furthermore, 
linear stability of the stationary and moving pulse solutions have been analyzed in 
P); see also EllEll- A bifurcation analysis of the self-sustained oscillations due to 
pulse dynamics near critical values of (j) can be found in Asymptotic analyses of 
the Gunn effect for the Kroemer model with a bistable g{u) can be found in |14l 1^. 

2.1. Bistable source g{u). The simplest case to consider is that of an N-shaped 
nonlinearity g{u) (the 'velocity-field characteristic') for u > 0, with a local maximum 
Qm = g{uM), um > 0, followed by a local minimum g^ = g{um) > 0, Um > um- Then 
J — g{u) may have up to three positive zeroes for J > 0, namely ui(J) < U2{J) < 
u^{J). For large enough L and K, gujuu < P < gm/um, and for in a certain 
subinterval of {uM,U3{gM)), there are stable time periodic solutions of l|2.1ll - 12.3(1 
of Gunn type. As shown in Fig. Il.lf a). while J{t) oscillates periodically, pulses of 
u{x,t) are created at a; = 0, move towards x = L, and disappear there. 

The analysis of the model is simple in the asymptotic limit 

(2.5) 0<e=y<l. 

In this limit H2.1|l - H2.2|l may be written as 

du du d^u J — n(u) 

\- K — — e ' 

ds dy dy"^ 

(2.7) / u{y,s)dy 



(2.6) ^ + ^^ = 4.,. 







where 



(2.8) 



y = ex, s = et. 
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Eq. (|2.6|l is a parabolic equation with fast reaction and slow diffusion terms. As 
e 0+, u(y, s) is typically a piecewise constant function taking on the order one 
values Ml (J) or M3(J) in intervals of length y = 0{1). The extrema of these intervals 
are typically moving internal layers. These layers are important because they bound 
pulses, and the self-sustained oscillation we want to describe is due to recycling and 
motion of pulses at the boundaries. A pulse is a region where u is U3(J), separated by 
moving wavefronts from two other regions where u is ui( J). There are two wavefronts 
bounding the pulse. In the backfront u increases from mi(J) to u^i J); this front moves 
with a speed c+(J). The forefront moves at speed C-{J), and in it u decreases from 
U3(J) to ui{J). Forefront and backfront are heteroclinic trajectories connecting the 
two saddles (wi(J), 0) and {u^{J), 0) in an appropriate phase plane (it, du/dx), where 
X = [y — Y{s)\/€ is a moving coordinate (x = at the wavefront and dY/ds — c±{J)). 
The instantaneous value of J{s) is determined by using the integral condition H2.7|l . 
Typically J obeys the simple equation 

(2.9) ^=A(J)Kc+(J)-n_c_(J)], 

where A[J) > is a known function of J, and n+ and n_ are the numbers of wave- 
fronts with increasing and decreasing u profiles, respectively. For high-field domains, 
n+ — 11- = 0, 1 [HIE]- The fixed points of Eq. H2.9|l correspond to the equal area rule 

(2.10) / [g{u) - c\du = Q 

J ui 

if n+ — ri-, or to possible plateaus in the shape of J{s) otherwise. Many questions on 
the stability of the pulses and their evolution can be simply answered by analyzing Eq. 
1)2. 9|) and using the asymptotic procedure explained in [51 . The description of pulse 
creation and annihilation at the boundaries requires a finer analysis, as explained in 

iniE]. 

2.2. Saturating source. If g{u) saturates, i.e. g{u) — >constant as u ^ oo, no 
u^iJ) exists and the previous construction is no longer possible. What is the correct 
asymptotic description of a pulse in this case? Let us anticipate the answer here. As 
Fig. Il.lf b) shows, a pulse is a traveling wave whose profile has a single maximum, 
and tends to u = ui[J) as a; — > ±cx). The pulse is bounded by a leading front 
(forefront, moving at speed c_) and a trailing front (backfront, moving at speed c+). 
Two parameters uniquely determine the pulse: its maximum height, Um, and the 
instantaneous value of J. Given Um and J, we can find the wavefronts enclosing 
the pulse by simple phase plane constructions. Consider the phase plane (u, du/dx), 
X — X — X{t), where X{t) is the instantaneous position of a wavefront and dX/dt 
its speed. There is a unique speed c+ = c+(J, u™) for which a separatrix issuing 
from the saddle (wi(J), 0) on the upper half plane reaches the u axis at {um, 0). This 
separatrix constitutes the backfront of the pulse, and a similar construction supplies 
its forefront moving at speed c_(J, «„). In general, c_|_ ^ c_, which implies that our 
pulse changes its size as it moves. How do we characterize the dynamics of pulses? 

Suppose that there is a single pulse moving in the sample. We need to determine 
the instantaneous values of J and u,„, for they completely characterize the pulse. 
The pulse width, I, changes as dl/dt = c_ — c^.. On the other hand, I may be 
determined by a line integral on the corresponding phase plane trajectories which 
form the pulse. Then I = 1{J, Um)- The dc bias condition yields a connection between 
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u„i and J, u,„ — U{J). Then the pulse width is a function of J only, ip{J) — 1{J, U{J)). 
Therefore since dl/ds = ip'(J) dJ/ds = c_ — c+, we get 

dJ _ c+{J,U{J))~c^{J,U{J)) 
^ ' ds " -^'(J) 

Typically the fixed point of this equation, J = J*, is such that c+ = c_ is a globally 
stable solution, so that J tends exponentially fast to J*. The corresponding pulse 
moves steadily without changing its size. This pulse is the homoclinic orbit in the 
phase plane, usually given by an equal-area rule and exhaustively studied by previous 
authors. Notice that the present construction explains why this steadily moving pulse 
is stable, thereby clarifying an old issue at the heart of the Gunn effect |2H1E1- When 
the pulse reaches the receiving contact, a different stage of the Gunn oscillation begins. 
This stage and others needed to fully describe the Gunn oscillation will be explained 
later. 

A subtle point is the following. Due to the integral condition l|2.2(l . the pulse 
height and width [in the (y, s) scales] are 0(e^2) ^ l and 0(e2) respectively, while 
outside the pulse, u — ui — 0(1) and J — 0(1). Thus our leading order approxima- 
tion for s) is not uniformly of the same order in space. Successive approximations 
of a single pulse lead to the following ansatz for u: 

u{y,s;e) ~ u^°\y,s;e) + e u'^^\y , s; e) . 

Here each u^j\y,s;e) may be of different order in e for different values of the space 
and time variables. However we shall impose that 

uniformly in y and s as e — > 0+. This situation results in a changing (self-adjusting) 
time scale for the evolution of J described by H2.11|l . See Ref. [7| for the description 
of a similar situation in combustion theory. 

3. Boundary layers and wavefronts. The model (12. 1|) - H2.4|l was introduced 
in order to argue that the asymptotics of the Gunn effect is universal within a class 
of models 6 . The nonlinearity g{u) was originally N-shaped: it had three branches 
for u > (with a maximum at um > and a minimum at Umin > um such that 
g{uM) > <?(umin) > 0, with g{u) — > oo as m ^ oo). In the present paper, we 
shall assume that g{u) is a smooth function having a single maximum at um > 0, 
g{uM) ~ gM > 0, such that g'{0) ^ (3 > 0, and lim.u^ooff(u) = a G (0, ^a/). To obtain 
explicit analytic expressions, we shall use a piecewise linear version of g{u), 

(3.1) g{u) — j3u9{uM — u) + a6{u — um) , 

where 9{x) is the Heaviside unit step function. See figure 13.11 where we have also 
shown a typical straight line J = u/ p which gives the value of u at the left boundary. 
Notice that the straight line intersects g{u) at {um , Pum), with a < puM < gM- We 
want to find stable solutions {u{x,t), J(t)) of Eqs. (|2.1|l and (|2.2|l under the boundary 
conditions H2.3|l and the initial condition (|2.4|) for appropriate positive values of the 
bias (j), in the asymptotic limit L — > oo. 
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Fig. 3.1. The nonlinearity function g{u) (solid), and the curve corresponding to the boundary 
conditions: J = u/p (dashed). The parameter values are <f) = 1.6, L = 200, K = 2.0, um = 2.0, 
/3 = 1.5, a = 1.0, p= 1.2. 



3.1. Outer limit and boundary layers. If < (f) < um, there is a single 
stationary solution of Eqs. H2.1|l - (|2.3|l which is easily constructed. Let y = x/ L = ex, 
< e <C 1. Introducing this scaling in Eq. H2.1|l yields the leading order equation 

(3.2) J-g{u) = Q, 

outside the boundary layers ai y — For < J < gM and < y < 1, (|3.2|) has the 
solutions ui{J) < U2{J). Of these we may choose u = ui(J), < ?; < 1, as the outer 
limit. Inserting it into H2.2|l . we find ui( J) + 0{e) — and then (|3.2|) yields 

(3.3) J^gW- 

For the piecewise linear g{u), J — (i<j), and ui{J) = J / [3. 

The boundary layers at j/ = 0, 1 are solutions of the problems 

(3.4) —^K—^J-g{U)^0, 0<e<oo, 

(3.5) U{Q)=pJ, U{oo) = ui{J). 

Here ^ = x = y / e for the injecting boundary layer aX y — Q, and = L — x = {l — y)/e 
for the receiving boundary layer at y = 1. The minus (resp. plus) sign in (|3.4|l 
corresponds to the injecting (resp. receiving) boundary. Clearly the shape of the 
unique solution of H3.4|) and (|3.5|l (for ^ — x) depends on whether J is smaller or larger 
than J = Jc = umI P- When < J < Jc, the boundary layer profile monotonically 
decreases from u = pj to ui = J/ [3. For the piecewise linear g{u), we have 



r 2 (pj - () 

(3.7) / (C/-Mi)de- ^ ^ 



^JW+APtK 
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However, ii J > J^, the boundary layer profile reaches a maximum before decreasing 
to ui — J/ (3. Numerical simulations show that in this case, the stationary solution 
of the model becomes unstable to Gunn-type oscillations. Given this occurs for 

(j) > Ul(Jc). 

3.2. Wavefronts. 

3.2.1. Bistable source. Let us first review how to compute the wavefronts 
when g{u) is N-shaped. Then H3.2|) has three solutions ui{J) < U2{J) < W3(J) for 
g{um) < J < g[uM)- As explained in Ref. |6j, the building blocks of the Gunn- 
oscillation asymptotics are wavefronts connecting ui{J) and M3(J). These wavefronts 
adjust themselves instantaneously to the value of J, as this unknown evolves on a 
slower time scale (see below). A wavefront centered at a; = X±{t) is a monotone 
function of x = — X± (t) such that 

m(— oo; c+) = ui( J), u(+oo; c+) = J) and 
m(— oo; c_) = J), u(+oo; c_) = iti( J). 

For the simple model used here, there is a relation between the wavefronts m(x; c+) 
and w(x;c_): 

Theorem 1. Let u(x;c±) be the wavefront satisfying 



(3.8) -(K-c±) — + J-g(u)^0, 

ax ax 

(3.9) u(— oo; c+) = ui(J), u(+oo; c+) = 1*3(7), 

(3.10) u(— 00; c_) = M3( J), u(+oo; c_) = ui(J), 

where x — x~X± (t) (X± is the position of the front at time t, determined by imposing 
u(0) — u'^ . dX±/dt ~ c±). Then we have 

(3.11) u(x;c_) = u(-x;c+), c+ + c^=2K. 



The proof is evident. This theorem shows that we only need to construct m(x; c+) 
and find c+ in order to have m(x; c_) and c_. 

3.2.2. Saturating source. Let now g{u) be a function with only two branches 
such as (|3.1|l . A wavefront is the only monotone trajectory connecting (ui(J),0) 
and a given point on the u axis (um,0). There is again a symmetry result for these 
wavefronts: 

Theorem 2. Let m(x;c+) be the wavefront satisfying iS.S\) . du/dx > for 
-00 < x = X - X+{t) < Xm, and 

u(— 00; c+) = ui( J), u{0;c+) = u^, 
du 

(3.12) u(Xm; C+) = Urn, Q^^^""' '^+) 

Here dX^/dt ~ c_|_, and 2xm = l{ J,Um) > is a function of J and Um- Then the 
wavefront satisfying J^S. du/dx < for — Xm < X = x~ X-{t) < 00, c_ = dX-/dt, 
and 

u{-Xrn;C-) ^ Um, ("Xm i C_ ) = 0, 

OX 

(3.13) w(0;c_) = u", m(+oo; c_) = Mi( J), 
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is such that holds. 

Again the proof is immediate. 

Let us now choose a certain Um = 
()2.2|l holds for a pulse made out of: 
• a backfront m(x; c+) at a; = 



U{J) for each (j)> um so that the bias condition 



x+{t), 



X+{t)\ and 



• a forefront u(— x; c+) at X- = X+ + 2xm- Now x = x — X-{t). 
The pulse u(x; c+) can be constructed explicitly for the piecewise linear g{u) of H3.1(l . 
If we choose — um, the solution of Eqs. I|3.8|l and l|3.9|) which is continuous and 
has a continuous first derivative at x = is 



(3.14) 
(3.15) 



m(x;c+) 
'«(x;c+) 



Ul{J) 



J — a 



X<0, 



UM 



K 



x- 



1 



X>0. 



Here mi(J) = J/f3, and 



(3.16) 



(3.17) 

(3.18) 
(3.19) 




if -c 

1 



K-c^ 



In 



(if-c+)2R 

A+(mm - ui)(-f'^ - c+) 



J - a 

Um = UM + — Xrn + 

K — c+ 



J — a 



If these expressions are inserted in the bias condition (|2.2|l . c+, Xm and Um may be 
determined as functions of J for a fixed (j). Figures IX^ and IX^ show the phase planes 
and the trajectories corresponding to u(x;c+) and u(x;c_), respectively, for a given 
value of J = 1.3. Fig. 14. II shows Xm and Um as functions of J. 

There are two important approximations to the wavefronts of Theorem 2, which 
yield either approximately triangular pulses (c+ 7^ c_, Um ^ 1) or the homoclinic 
pulse (c+ = c_ = K). These two limiting cases are described in detail in Appendix 

HI 

4. Pulse Dynamics. 

4.1. One pulse far from the boundaries. Let us consider a single pulse 
moving far from the boundaries, as described in the previous Section. Its height 
and width are established by imposing the integral condition H2.2|l . The result will 
show that the pulse is a tall and narrow moving layer which changes size as it moves. 
At certain stages of the periodic Gunn oscillation (see below), the pulse height is 
u — 0(e^2), its width is Ay = 0(e2), so that the pulse excess area (in y space 
units) is 0(1). The inner core of the pulse contributes an order 1 amount to the bias 
whereas its exponential tails approaching ui(J) yield an 0(e) correction to its excess 
area. Thus we suggest a decomposition of the solution u into an outer solution valid 
outside the pulse and boundary layers, and inner solutions comprising the front (s) and 
boundary layers. In this Section, we shall calculate the leading order term in each 
asymptotic expansion explicitly. The first correction to these results can be found in 
Appendix ^ 
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Fig. 3.2. Phase plane (u, du/dx) and trajectory corresponding to n(x;c+) for J = 1.3, c+ = 
1.8359. Other parameter values as in Figure rOl 




Fig. 3.3. Phase plane (u, du/dx) '^'"■'^ trajectory corresponding to n(x;c_) for J = 1.3, c_ = 
2.1641. Other parameter values as in Fiaure T3.1\ 



4.1.1. Outer solution. The outer solution is 

(4.1) = u(0)(y, s) + eu^^\v, s) + 0(e2) 

(4.2) J = J("'(s)+ej(i'(s) + 0(e2). 

Here = 0(1) yields an order 1 contribution to the integral condition H2.2|l . of the 
same order as that provided by integration of the excess area of the pulse inner core, 
eu'^' yields an order e contribution to H2.2|l . of the same order as that provided by 
integration of the excess area of the pulse tails. Inserting this Ansatz into (|2.t)|) . we 
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obtain 

(4.3) J("' = 0, 

(4.4) J«-,V°V^^^ = ^+i^^. 
Solving and llOl yields 

7(0) 

(4.5) jo)^^^(j(o))^:^^ 

7(1) 1 d./'"' , , /nxs 

J^°) and J^^) will be found later from the integral constraint. 

4.1.2. Two-term description of the pulse. The pulse described in the pre- 
vious Section may be considered a moving inner layer solution. Its height is much 
larger than 1 and its width much smaller than 1 (the precise orders will be determined 
later). We shall assume 

(4.7) u™"- ~ (y, s;e) + e P^'^ (y, s; e), 
where P*^"^ is the pulse solution of H3.8|) described in Theorem 2: 

P(°) (y, s; e) = u{x - X+; c+) 9[xm - {x - X+)] 

(4.8) +u{X+ + 2xm - X- c+) e[x -X+- Xm]. 

P^"-* {y, s; e) and P*^^^ (y, s; e) depend on e in such a way that 
P(i)(y,s;e) 



P(0)(y,s;e) 



= 0(1), as 



uniformly in y, s. The inner expansion H4.7|l is chosen so that the pulse inner core 
in P^°) yields an order 1 contribution to the bias condition whereas its tails together 
with the inner core of eP^^-* yield an 0(e) contribution. The latter is of the same 
order as the contribution of the outer solution, eu^^-*, to the bias. 
In Eq. H4.8|l . the location of the fronts and their velocities are: 

(t) + rf«(t), 

(4.9) c+^cf+ec^^\ 

and similarly for X_ and c_. Inserting (|4.2|l and H4.7|l - (|4.9|l in p.lfl and (|3.9|l . we 
obtain for the trailing front 

(4.10) -{K- cf) ^ + J(°) - .g(p(°)) = , 
OX OX 

p(°)(-oo;cf ) = = u,, P(")(0;cf ) = um. 

(4.11) P(°Hxm;c?^)-«m, ^(x™;cf) = 0. 
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whose solution is 14. 8|) . and 

(4.13) pW(-oo;cr) = ^(x™;cV"^) = 0. 

The leading front of the pulse obeys similar expressions. The correction P^'^^ is cal- 
culated in Appendix^ in which explicit formulas for piecewise linear g{u) are given. 

4.2. General equation for Using Theorem 2, we can calculate the bias 

condition for a single pulse moving far from the boundaries as 

Jo 

1) ^.(0) 
2 



J-oo 9l 



(4.14) 



+e [UL{0-ui]d^ 
Jo 

' r[UR{0-u,]dC + O{e'). 
Jo 



Here the bias is the sum of the areas of the regions inside and outside the moving 
pulse. The leading order contributions to these areas are the first two terms on the 
right side, which are of order 1. They are: (i) the leading order contribution of the 
outer solution and (it) the leading order contribution of the inner core of the pulse. 
The other terms are 0(e) and correspond to: 

• (i) the tails of the pulse to leading order, 



2e 



oo 



• (ii) the second order contribution to the outer solution, 

• (iii) the second order contribution to the area of the inner core of the pulse, 
2e^ /o^"" P^^^dx = 0{e) (we have ignored a much smaller term of order f^Xm)^ 

• (iv) the layer at the left boundary, Eqs. H3.6|l and l|3.7l) : 



{Ul - ui) dS, = — — , 

^K^ +A(3- K 



The area of the injecting (left) boundary layer becomes of order 1 when a new pulse 
is being shed; otherwise it is of order e, as indicated above. 
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If the pulse is far from the boundaries only the first two terms are of order 1, and 
we have 



2(if-cf ) 



K — c 



(0) 



0(6) 



Ul + 2e 



Xm{uM - Ul - -B+) + 



2(i^-cV°^) 



(4.15) 



B+\+{uM - Ul) 



J(") - a 



0{e). 



Now we can proceed as sketched in Section |21 Eq. H3.18|l allows us to obtain c^'^^ 
as a function of J^"^ and Xm'- 



(4.16) 



for a fixed value of </>. Inserting this function in (|4.15|) . we can determine Xm as a 
function of Then yields function of J*^"^. The results are certain 

functions: 



(4.17) 



5(J(0)) 



that have been depicted in Fig. 14.11 
Eq. lfTTT|l : 



Time differentiation of 2xm = fiJ^'^'^) yields 



(4.18) 



(J(0),t/(J(0)))-X 



-(/5'(J(0)) 



which describes the time evolution of J^^\ Provided that J^°)(0) > J*, J*-"^ decreases 
exponentially fast to J* such that c^|"'' — cf"'' = if. The pulse then moves at constant 
— J* and speed K, and it is a homoclinic orbit of the phase plane 1)3. 8|l with 
c = if, as shown in Fig. 14.21 

Figure 14.31 compares the leading order asymptotic solution with the numerical 
solution. The upper part shows the time evolution of J{t) and the lower part u{x, t) at 
the times marked in the upper figure. Notice that initially J{t) decreases exponentially 
fast to J*. Consider an instant (point 6) at which the wave is fully developed and 
far from the boundaries. We observe that the profile of the asymptotic solution 
has a larger height and is slightly thinner than the numerical solution. Why? We 
have neglected 0(e) terms (boundary layers, wave tails, . . . ) when calculating the 
integral condition with the leading order asymptotic solution. Then J^^^ (t) is slightly 
larger than the numerically calculated J{t). The asymptotic profile fully agrees with 
the homoclinic solution described in Appendix IXI bv l|A.10l) . (|A.lip and (|A.12p (also 
calculated excluding order e effects). 

The previous ideas are correct if we can show that evolves on a slow time 
scale, say a — y/e{t — to) or t — e{t — to). To see this we should analyze ()4.15(l and the 
previous equations with a little more care. We shall show that there are two limiting 
cases for which the pulse can be easily calculated and H2.11|l explicitly obtained to 
leading order: 
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15.0 




Fig. 4.1. Half-width Xm and maximum height Um of the single pulse as a function of J for 
the same parameter i^nht.ps as in Fimirps \S 1\ tn \s .SI 




Fig. 4.2. Phase plane (u, du/dx) <^^d, homoclinic orbit for J* = 1.13182, c = 2. Other 
parameter values are as in FiaureUl\ 



• - J* = 0(1), K - c^°^ = 0(1) and - a)xm > 1- The limiting pulse 
is a triangular wave formed by pieces of phase plane trajectories which do not 
tend exponentially to infinity as x — * ±oo. 

• [K — c^'') ^ (J*^"^ — a) ^ 1, Xm ^ 1. The limiting pulse is the homoclinic 
trajectory with c — K. 

In the first case, J*^"^ evolves on the slow time scale a = \f(.[t — to). In the second 
case, the time over which J^*^) varies appreciably is o(e~^). One period of the Gunn 
oscillation contains stages during which these limiting cases are good approximations 
to H2.11|l . In fact the Gunn oscillation may be considered as transitions from one 
limiting case to the other depending on the number of pulses existing at the time. 
Then changes slowly with time which justifies our wavefront constructions. 

Appendix [CI analyzes the dynamics of pulses corresponding to the two limiting 
cases of triangular and homoclinic pulses. At this point, we may call attention to two 




Fig. 4.3. Upper: time evolution of J. Times marked correspond to; (1) wave reaches receiving 
contact, (2) J reaches Jc (3) Xm < XL (A) wave completely exits (5) J < Jc (6) fully developed 
wave. Lower: profiles of u(x,t) at the times marked in the upper figure. 



peculiar features of our results: 

1. Our leading-order approximation for the solution is of order one outside a 
pulse, but it is much larger inside it: of order e~2 . 

2. The proper time scale for the variation of J is i if the pulse may be approxi- 
mated by a triangular wave. It is slower as the pulse approaches a homoclinic 
pulse. Both these approximations are limits of the same equation, (|2.11|l . 
Thus this equation contains more than one asymptotic limit and its time 
scale is changing as time changes. 

5. Pulse dynamics near the boundaries. In the previous Section, we de- 
scribed the dynamics of a single pulse far from the boundaries. Basically the pulse 
approaches the homoclinic pulse on a slow time scale: at each time on the scale t 
the pulse follows adiabatically the instantaneous value of J. J changes on the scale 
(T = ^/et or even more slowly as described by Eq. (|2.11|) . where (/j is a function of J 
and e {^p{J) is of order e~2 or larger). In this Section we shall describe what occurs 
when the pulse reaches the receiving boundary at a; = L and beyond, until a period 
of the Gunn oscillation is completed. 

5.1. Pulse disappearing through the receiving boundary. Let us assume 
that a single pulse has reached its asymptotic shape (J*^"-* — J*) and it advances 
with speed K until its forefront reaches X- = L at time ti. Afterwards, it begins 
leaving the sample. As time elapses the wave exits through the receiving boundary 
and therefore the area under the wave decreases. Since the total area has to satisfy 
the bias condition, this loss of area has to be compensated by a corresponding increase 
in ui so that equation H2.2|l still holds. 

Let us denote by ul the value of the pulse inner solution P*-*^) (y, s) at the receiving 
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boundary, ul is obtained from equation H3.15|l for u{x ~ xi"-';c'^-') — u{x'^^ — 
x; c^"* ) = u{2xm ~ X\ c^"* ) when x ~ 2xm = L — xl^' . The corresponding argument of 
- x\ c^"^) is XL = - L > 0. When < xl < Xm, the bias condition ijT^ 

yields 

(5.1) = + $W(J,cf +e$«(J,cf ,x™,XL, + 0{e^), 

<i>(0)(J,cV\x™,XL) , , (0). 1 , . (0). , , 



(5.2) 



(2Xm - xl) (um -Ul- B+) + /'^ fx™ - 4 



$W(J,cf,x™,XL,/^^)= f {P^"-^-u,)dx+ ^^''-,f' 



oo 9l 



2 



/■X™ (-Xt /-oo 

(5.3) +2e / p(^)dx-e P^^Ux + e [Ul (C) - «i] d^, 

Jo Jo Jo 

instead of 14.14(1 . In this equation and ^^^^ are of order 1 because the integrations 
of P*^°) and p(^) over the inner core of the pulse are of order e""'^. Eq. H5.1|) yields 

(5.4) = 7.i + $(")(j(°),cf ,x™,Xl), 

(5.5) a>«(j(°),c^"\xL,x™,J<^') = o. 

We shall now find the evolution equation for J*-"-* by a procedure similar to that 
used to find H4.18|l . The right side of ((5.4|l depends on J^^\ c-^\ Xm and xl- Xm is a 
function of j'"' and c^^ given by H3.18|l . As wavefront velocity c"^\ we shall use the 
function of J^°) (for a fixed value of </>) that was determined at the end of Section O 
Then the right side of 1)5. 4|l is a function of J^*^^ and xl only (for fixed 0): 

</) = S(J(°),xl) = wi(J^"^) 

(5.6) +<i>(°)(j(°),cV"V^°'),X.^(j(°\c?)(j("))),XL). 

Xl can be explicitly calculated from 



(5.7) ^^cL°^=2i^-cf, XL(ti) = 0, 



where ' is our known function of J^^'^ . We can obtain a closed system of equations 
for XL and by differentiating ((5.6(1 with respect to time and then using 1(5.7(1 . The 
result is 



<^-«' 



Here we have used that (15.2(1 and ((5.6(1 imply dB/dxL = —{ul — ui). 

The time evolution of j'"^ is found by solving this equation while the wave dis- 
appears through the receiving contact. Having found the solution to leading order, 
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(|5.5|) yields the correction j'^^K Numerical solution of H5.1|) to H5.8|l shows that 
increases with time, as shown in the region between times 1 and 2 in figure 1^751 Notice 
that this increase agrees with the numerical solution of H2.1I2.4|I . 

Depending on the bias 0, one of the following two events may happen first: 

• (i) J*-*^-* reaches Jc, or 

• (ii) XL = Xm- 

In both cases the stage described by the previous equations ends. In case (i), a new 
wave is created at a; = 0, whereas in case (ii) (15.111 should be changed to 

(5.9) ^ = «i + *(°)(j(°),cf ,XL) + e*W(j("),cV°\xL,jW)+0(£'), 

( Jf") , cV" ' ,XL) = e Hx; c^°^ ) - ^.i] dx 

{j("^-a)xl 







(5.10) 



(0) rXL{uM-Ui~B+) + — 



2{K^c';') K~c 
*W(J,cf ,XL,J«)= r Hx;cf)-u,]dx+ '^''^-,f^' 

J-oo 9l 

/•XL i-OO 

(5.11) +e / P^^Ux + e [UL{£)-ui]di 



(0) 



JO JO 

where now XL = L — X,^^-* > and 

(5.12) XL{t2)^Xm{t2). 

Here ^2 is the time at which the maximum of the pulse reaches x — L (equivalently 
XL — Xrn in the previous stage). Similar arguments to those used in the previous 
stage lead to, 



/i 

(5.13) I ^+9JW ) ^'^""^ 



dt 



-+ 



where ¥^\J^''\xl) = «'(°)(^^°^4"^(^^°^),Xl)- 

This stage lasts until reaches Jc and a new pulse is shed. If the bias is 
small enough, the solution of these equations may be such that J*^°^ never reaches 
Jc and it eventually decreases with time. In such case, the pulse exits and leaves a 
stable stationary state in its wake after xl = 0. Notice that setting xl = Xm = 
in H3.18|l implies a front velocity ljA.7|l . Thus the velocity of the disappearing wave- 
front approaches that of the triangular wave, although the shapes of the respective 
wavefronts may differ greatly. As in the previous stage, we find J*^^^ by solving the 
equation ^^'i {J^'^\cf\xL, J^^^) = 0. 

5.2. Pulse shedding at the injecting boundary. If J{t) reaches Jc at t — t2, 
the boundary layer becomes unstable and a new pulse starts being shed at the in- 
jecting boundary. The boundary layer profile, U{x, t) solves the following semiinfinite 
integrodifferential problem: 

(5.15) U{0,a)^pJ, 
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whose solution exhibits an explosion-type behavior and a rapid growth of the area 
enclosed by the boundary layer which can no longer be neglected. This increase in 
area is to be compared with the area released by the disappearing pulse at x = L. 
Initially, the area released is larger and this net area loss has to be compensated by 
an equal increase in the area of the outer solution, so that J continues to increase 
although at a slower rate. After a short time, the growth of the boundary layer is 
larger than the area released, so that J reaches a maximum and starts decreasing. 
The structure of the injecting boundary layer when J > Jc is as follows: 

1. u = u{x; J, Urn) is quasi-stationary for < x < Xm such that H3.4|l holds with 
u{0;J,Uni) = pJ, u{Xm; J,Um) = Urn, with du{Xm; J,Um)/dx = 0. The 
boundary layer solution is the trajectory of the phase plane corresponding to 
(|3.4|l which leaves the vertical line u = pJ a.t x = and intersects the u axis 
at u = Um- Notice that this trajectory is uniquely determined by giving J 
and Um- 

2. For X > Xm, the boundary layer is a wavefront of type w(x;c_) moving at 
speed C- . This speed and the forefront are uniquely determined by J and 
Urn- u = C+), ^ = — X, with C+ = 2K — m(— oo; C+) = ui{J), 
u(0; C+) = UM, C+) = Um, du{U; C+)/d^ = 0. 

The bias condition (including the injecting boundary layer) is 

(5.16) = + (J, Xrn, XL)+e ^'"^ (t) + 0(e), 

(5.17) = /"[C/(e)-C/oo]dC, 

provided that the maximum of the old pulse has leaved the sample. Time differenti- 
ation of (|5.16|l yields 

^'■''^ [p + ojm)^- + ^ - ^V^- 

Equations and lfOH|l . together with and l|lT7|l . constitute 

a closed system of equations for the unknowns J, xl, and Xm- During this stage, J 
initially increases and then decreases until it either reaches again Jc, or the old wave 
completely disappears. In the latter case, the evolution of J{t) is still given by equation 
H5.18|l without the last term which represents the area lost by the disappearing wave. 
The evolution during this stage can be observed between points 3 and 4 in Figure IT!^ 
Also, Figure I5T1 shows in detail the evolution of the boundary layer profile from 
the time that J > Jc- Initially, u grows at the injecting boundary because J increases. 
Furthermore, the slope du{0, t)/dx increases and it becomes positive at a certain time. 
Then a wave-like structure is created. The leading front of this wave moves away from 
the boundary while the back- front is attached to it. The wave continues its growth 
until the time when J becomes again smaller than J^,. Then, the slope at the boundary 
becomes negative, the wave dettaches and moves away from the boundary as a solitary 
wave. At that time, the quasistationary part of the boundary layer which joins x = 
to the maximum at x — Xm becomes the backfront of a detached pulse. Then, we 
again have the same equation (|2.11|l describing the first stage, and a cycle of the Gunn 
oscillation has been completed. The time evolution during this stage can be observed 
in Figure ^31 for times larger than ts. 

6. Conclusions. In this paper we have asymptotically described one period of 
a Gunn-type oscillation in a simple model. The nonlinearity of the model is such 
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Fig. 5.1. Boundary layer profile evolution during the shedding stage. 

that at most two constant solutions are possible for each value of J (the current- 
like unknown). The model consists of a parabolic equation for the field-like unknown, 
u{x, t), and an integral constraint (bias condition) which determines J{t). Appropriate 
boundary and initial conditions are imposed. The key new idea of our analysis is that 
a pulse that changes shape as it advances may be constructed by fixing only two 
parameters: J and the pulse maximum, Um. If the pulse width is small compared 
to the sample length, L, J and Um change on a slow time scale. The trailing front 
of the pulse is part of a separatrix joining a saddle point to (u™, 0) with du/dx > 
on the {u, du/dx) phase plane. The initial and final points determine the backfront 
speed c+ as a function of J and Um- Similarly the forefront of a pulse is constructed 
and its speed c_ determined. Then equations for J and obtained by time- 

differentiation of the bias condition and of the pulse width. The time derivative of the 
later is c_ — c+ and J tends towards a fixed value corresponding to a rigidly moving 
pulse with c_ = c+. Other stages of a Gunn oscillation including wave creation and 
annihilation at the boundaries are analyzed by similar methods. Our theory compares 
well with direct numerical simulations. 

We have found an asymptotic theory of the 'Gunn effect' in a simple piecewise 
linear model, whose main step is a construction of pulses and a derivation of an 
equation for the 'current'. An analysis of the stability of these solutions is an open 
problem, although there is some work on this problem in the related Kroemer model 
[3 125j . That the profile of pulses can become oscillatory for appropriate parameter 
values has been shown by Sun et al for some reaction-diffusion models We expect 
that the present method yield results independent of the model equations within a 
class thereof displaying the Gunn instability [BJ |^ \'2'2i 1251 I26j . Studies of other 
systems that can be understood by the dynamics of pulses are in progress. 

Appendix A. Limiting cases. 

A.l. Triangular pulses. The bias condition H2.2|l determines the orders of mag- 
nitude of Xrn and Um in terms of the small parameter e. Let us assume that {K — c+) 
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and (J - a) are 0(1), whereas Xm > 1- Then iflTTTjl . (|TTH|) and ifTT^ imply that 

J — a X+{um — ui) 

Urn = UM + -r^ Xm 

A — C+ A — C+ 

/ A 1 N J - a 

A.l Xm. 

The wavefront u{x: c+) is given by H3.14|l for x < and 



J — a 

(A.2) u^U.^ + - {x-Xrn) + B+ 

K — c+ 

/ . „N J ~ a 
A.3 ~- X, 

A — C-|_ 

for X > Oj where IjA.lfl has been used. The bias condition (|2.2|) then yields 
(A.4) ^-^XL 



(A.5) Xm 



- Ml) {K ~ c+) 
e(J- a) 



These equations imply that Xm and 0(e 2 )^ while the proper time scale over 

which J varies is 

We can obtain (jA.ip - ljA.5|) directly from the Eq. (|3.8|l and the bias condition. If 
Xm 2> 1, so is Um- Then we shall select uniquely the wavefront solution of Eq. (j3.8|l . 
u(x; c+), so that it satisfies u{—oo\ c+) = ui(J) and it does not grow exponentially as 
X +0O. Similarly, u(x;c_) = u(— x;c+) does not grow exponentially as x — * — oo, 
and it satisfies u(+cx3; c_) = ui{J). As in Theorem 1, we still have c+ + c_ = 2K. 

u(x\c+) satisfies Eq. (|3.14|l and 

(A.6) u(x; c+) = UM + ^7 — — X, if X > 0. 

A — C+ 

Continuity of du/ dx at x = directly yields 



J\ , J-a 



(3 J ^ X - c+ 



which in turn implies 



(A.7) c+ = K ^ " > 0. 



\j P (uM - Ui) (^UM - f ) 



Figure lA.ll compares the actual values of c± with the approximation ljA.7|) . Notice 
that both lines are resonably close for J sufficiently higher than J*. 

We can now form a pulse by joining the backfront u(x; c+), x — ^ ^ X+ < Xm 
to the forefront u(2xm — x;c+). This pulse asymptotically approaches an isosceles 
triangle of basis {X- — A"+) — 2xm and height approximately given by 



Um = ""(Xm; C+) 

(A.8) 



2 

(J-a) (a:_ - X+) 
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Fig. A.l. Velocities of the backfront, c^, and of the forefront, c_, as functions of J for (f) = 1.6 
and the same parameter values as in Fiaure W^ We have also shown the corresponding approximate 
values c± obtained for the triangular pulse described in the text. 




Fig. A. 2. Shape of the pulse for J = 1.3 and the parameter values of Figure Ty. 11 Also shown 
are the corresponding triangular and homoclinic pulses. 



See Fig. rOl which compares the triangular pulse to the real pulse and the homoclinic 
pulse for the same values of J and cj). Here we have used (IA.6|) and assumed that 
Xm (the location of the maximum, equal to the pulse half-width) is very large. To be 
precise, we assume that {K — c_|_) = 0(1) and that (J — a) Xm ^ 1- 

Notice that the way we have constructed m(x; c+) is immediately applicable to a 
general smooth nonlinearity g(u) of the same type. We have thus the following result: 
Result. Let the characteristic curve g{u) be an odd function of u with a positive local 
maximum after which it monotonically decays to a positive constant as u +oo. 
The approximate backfront u(x;c+) is the unique solution of Eq. I^!^.<^} which, for an 
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appropriate choice of the velocity c+, satisfies u(— oo;c+) = Wi(J) and it does not 
grow exponentially as x ^ +00. 

A. 2. The homoclinic pulse. If c+ = c_ = K, the pulse is a homoclinic orbit 
of the phase plane H3.8|l : 



(A.9) 



— + J-g(u) = 0, 



with u(±oo) = ui(J) = J//3. See Fig.E3 Here C = x-Xo, and Xq = {X+ + X_)/2 
is the location of the maximum of the pulse, Um. The solution is 



J 



+ [ UM - 



(A.IO) 
(All) 

(A.12) 



.-v^(ICI-Co) 



for Id > Co, 



m(C) = "m - ^^-^ for ICI < Co, 



Um = Um + 



J — a 



Co is determined by imposing continuity of du/d(^ at C = iCo- 

VJ3 (um - 



(A.13) 



Co 



J — a 

Now we may find J — J* from the bias condition. The result is 



(A.14) 



(A.15) 



(A.16) 



J* =a + /33 



Co = r' 



\ 2e [um - f 



'A 

2e 



0(e), 



o(i). 



0(1). 



Fig. IA.3I compares J* to the approximation (|A.14|) . Notice that a simple phase plane 
argument indicates that J* obeys the following equal-area rule: 



(A.17) 



r 



g{u) du , 



where m,„ is given by Eq. (|XT6^ . Fig.lXl compares the actual pulse, the homoclinic 
pulse with c — K and the triangular wave for the same values of J and 4>. 

Appendix B. Explicit calculation of F^^^ and c*^^\ First of aU, P^^^ = 

is a solution of the homogeneous problem H4.12|) (with zero right 
hand side) and H4.13|l . The solvability conditions for the nonhomogeneous problem 
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Fig. a. 3. Value J = J* for the homoclinic pulse as a function of the bias rf>. The dashed curve 
corresponds to the approximation \A.14f - 



yield 

(B.l) 
(B.2) 

Let us define 
(B.3) 



p(0) p(l) „ p(0) p(l) 
X X XX 



p(fi) p{Q) _ j(l) _ c^l)p^O) 



(0), ap(«) 



dx 



Inserting (|R3)| in (g^S, 

we obtain an equation which can be solved by two quadra- 
tures. We find 



(B.4) ^ = 



dx 



(0)2 



(K-cf )xp(0) \p{0) _ j(l) _ ^(1) p(0) 



dx- 



Here subscripts indicate partial derivatives with respect to the corresponding variable. 
It is easy to check that this expression satisfies HB.1|) . Further integration yields (for 
X<0) ^ 



P'^X+{uM - ui) 



(B.5^- 



t(0) 



A 



+s 



+ c 



(1) 
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where g is a constant. For x > 0, P^x = — c^|"'')^e*-^ ^+ ■'-^ and we can simplify 

the expression for p by integrating by parts. The result is 



1 



P = 



p(0) 
X 



-(i^-cf )x 



-(K-c<"))x 



(B.6) [pi") - J(i) - c«Pf )] rfx} + 

Continuity of p at x = yields 



(B.7) 



q = 



pig 



X+f3{uM - Ml) 



where Iq is the following integral: 



g-(K-cf ) X p(0) rp(0) j(l) _ c(l)p(0) 



whose value can be computed as, 



(B.8) 
where 



h = h 



00 



he c 



(1) 



^00 

h.J 

he 
3(1) 



(0) 



{UM - Ml) C+s 



(0)^2 



4/3 



Continuity of P;^ at x = implies 



(B.9) 



(0) K(0+)-p;,(0-)] = 

-p(o) fp("'(o+)-p(")(o- 



The second argument has been omitted in all the functions in this formula. The 
jump discontinuity of the second derivative ^-t x = implies that has also a 



jump discontinuity at x = 0. Substituting (|B.5p and (jB.6p in HB.9|I . we obtain 



(B.IO) 



QB+{K ~ cff = q\\{uM - ui) - 



J^°^ , ( \\ (1) , Ks{uM - Ul) 

+ [UM - ui)\+c\' J 



1 
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We can obtain Q from ljB.7|l . (|B.8|I and ljB.10|) . After some algebra, the result is 



(B.ll) 
with 



(1) 



+ 



AO) 

J s 

T 



+ ^00 



X+ [UM - Ui) c^l\ 



+ 



B+{K-cf Y [(if ~cV°^)2+4/?]3/2 



vQc = 



\+{uM-Ul) 



= A^(uM-ui) - B+{K- 



The velocity correction c^-* can be obtained from the condition 9P(^'/9x = at 
X = Xm- Thus 



Q 



B+{K-c 



{0)^ 



(B.12) 



Inserting IjB.llj) in this equation, we get c^^ . After simplification, the result is 



(B.13) 
with 



c+o + C+J 




Z2 



(3{um - ui) 



+ (j(") - a) zi + 



C+J 



-B+Qo{K-c 
-Z2- B+Qj{K ~c 



(0)n2 



(0)n2 



J(0)_ 



7 = 



is: - c 



(0) 



Z2 - P+(A- - )x™ + P+Qc(i^ - cV°^)2 



Zl 



J(0) - a 



_ exp [{K - c)_')xm]- I A+ (mm - mi ) 

^2 = — ^ 



(if-c?)) 



J(0) 



a 
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Here we have used (|3.16|) - (|3.18|l to simplify the result. J^^^ is found from the 
equation: 



= 2 
+2e 



9i 

{UL{0-ui)d^ 



(B.14) 





oo 



(C/fl(0-tii)de, 



after substitution of (|4.14|) . H4.15|l and (|B.13|I . 

Appendix C. Pulse Dynamics: Limiting cases. 

C.l. Triangular pulse. If the limiting pulse is triangular, the approximate 
evolution equation for J may be obtained by time-differentiating IjA.Sp and using 
iTOIl : 

(C.l) 
(C.2) 



B{J) 



^ = -^eBiJ) [K-c+{J)], 
4(3{Pum - a)i{(j) - ui)i {um - ui)i 



2um ^ ^ 

For typical values of the parameters such as those in Fig. 13.11 B{J) > 0, so that J 
tends to the solution of c+(J) — K. Triangular pulses are good approximations for 
J sufficiently large, which means that the solution of IjC.ip decreases towards J — a 
according to (|A.7|I . Of course before this value can be reached, (|C.1|I ceases to be 
valid and we revert to the general equation for J, H2.11|l . whose fixed point is J*. 

C.2. Homoclinic pulse. Let us now assume that {K — c+) <C ( J — a) <C 1, 
whereas Xm > 1- Then (|3.17|l and l|3.18|l imply that 



(C.3) 



J 



2(J-a)2 



Here A+ ~ ^/]3. Notice that ljC.3|) becomes (|A.13|I as {K — c+) ^ 0. We now insert 
this approximation in H3.19|l after H3.17|l has been substituted. The result is 



(C.4) 



1 



2{J-a) 

The bias condition (|4.15|) may now be approximated by using (|C.3|I and 
obtain 



1) to 



(j)-ui 2(32[uM-uif 



Then 
(C.5) 



(J -a) - 



' 2e(uM - uif 
3(0 - ui) 



2e (uM - 
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Inserting HC.5|I in HC.3|I and HC.4|I . we obtain 

a-i / 3(0 -til) 



(C.6) 



2e {uM - ui) ' 



(C.7) ^ V ^ '•'^ ^ ^ ■ 

Equations ljC.5|) - IjC.Tp are the same as (IA.14|) - (|A.16|) for the homochnic pulse. 

As explained before, dxm/dt = (c_ — c+)/2 = K — c+. Therefore the derivative 
of l)C.6|l with respect to time yield 

.^o^ dJ [8^ I3i {(f)-ui)iiuM -ui)i 

(C.8j — \— [K~c+). 

at V o Um — <P 

This equation has the same form as (|2.11|l , and it shows that the unknown J{t) varies 
on a slow time scale t = 0{l/y^e{K — c+). In the present hmit, (X — c+) <^ ( J — a) = 
0{y/e), so that the corresponding time scale is slower than t ^ et. 



A glance to Eq. (|C.8p shows that J decreases exponentially fast to J* such that 
c+ = c_ = K. The resulting pulse is the homoclinic orbit of the phase plane ()3.8|l 
with c = K described in the previous Section. 
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